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Quantum filter for a class of non-Markovian 

quantum systems 
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Abstract —In this paper we present a Markovian representation 
approach to constructing quantum filters for a class of non- 
Markovian quantum systems disturbed by Lorentzian noise. 
An ancillary system is introduced to convert white noise into 
Lorentzian noise which is injected into a principal system via 
a direct interaction. The resulting dynamics of the principal 
system are non-Markovian, which are driven by the Lorentzian 
noise. By probing the principal system, a quantum filter for the 
augmented system can be derived from standard theory, where 
the conditional state of the principal system can be obtained 
by tracing out the ancillary system. An example is provided to 
illustrate the non-Markovian dynamics of the principal system. 

I. Introduction 

The rapid growth of quantum information technology 
greatly boosts the development of quantum control theory Cl 
which in turn makes such technology more reliable. For exam¬ 
ple, quantum filtering IT] and quantum feedback control |2l, 
0, 0 have been applied in stabilizing quantum information 
carriers Q, in systematically designing photonic computation 
circuits 0, and in enhancing the performance of quantum 
metrology Q. 

The quantum plants that most existing works are concerned 
with are Markovian, and are disturbed by white noise only fS], 
E). However, many problems of interest involve complicated 
environmental influences, e.g., colored noise cni, C2, which 
if not taken into account may lead to degraded performance 
of estimation and control schemes. Similar modelling issues 
arise for classical systems IfT^ . and a common approach in 
control engineering to account for the non-Markovian effects 
of colored noise is to augment the system with a whitening 
filter 03. 

In this paper we consider non-Markovian quantum systems 
by representing them in a larger Markovian system framework. 
If we assume that the principal system of interest is defined 
on a Hilbert space 1), and the noise on a Fock space we 
introduce an ancillary system defined on the Hilbert space 
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()o converting white noise into colored noise to model the 
internal modes of the non-Markovian environment, whose 
structure determines the spectrum of the colored noise. Then 
the augmented system has a Markovian evolution on the 
augmented Hilbert space f) (8) f)o ® 5^- Such an approach 
was proposed in IH to model a non-Markovian system 
and was named as a pseudo-mode method later Ifl5l . This 
approach is also applied to model energy transfer process 
in photosynthetic complexes na. Similarly, the augmented 
system model can be realized by a quantum collision model, 
where the spectrum of the noise is implicitly determined by 
the ancillary system ini. Also, a hierarchy equation approach 
has been adopted to describe the dynamics of non-Markovian 
quantum systems na, where parts of the equations describe 
the pseudo-mode dynamics. This has been applied to the 
indirect continuous measurement of a non-Markovian quantum 
system m. However, this pseudo-mode approach has not 
been systematically described so as to be compatible with 
quantum control theory, e.g., quantum hltering theory. 

For concreteness, in this paper, we assume that the noise 
has a Lorentzian spectrum and employ a quantum stochastic 
differential equation (QSDE) approach to represent the col¬ 
ored noise via an internal mode of the environment. Due to 
the prevalence of Lorentzian noise in some solid-state sys¬ 
tems M, El, it is useful to model Lorentzian noise in quan¬ 
tum control applications. This colored noise is injected into 
the principal system of interest via a direct interaction such 
that the dynamics of the principal system can be described by 
a quantum stochastic integral differential equation (QSIDE), 
i.e., a non-Markovian Langevin equation. In addition, the 
augmented model of the non-Markovian quantum system can 
be conveniently described by an {S,L,H) description in an 
extended Hilbert space which is compatible with quantum 
filtering theory. The quantum filter for the non-Markovian 
quantum system can be constructed by injecting a probing field 
into the principal system. Due to the output field satisfying 
a non-demolition condition, the augmented system state can 
be estimated by the filter, with which the non-Markovian 
dynamics of the principal system can be obtained by tracing 
out the ancillary system. We note however that more general 
noise spectra can be considered in an analogous manner (more 
details will be presented in a future paper). 

The paper is organized as follows. In Section [III based on 
a Markovian quantum system model, we introduce a colored 
noise model to model Lorentzian noise. In Section |III1 we 
show that the principal system disturbed by the noise from 
the noise model satisfies a QSIDE. In Section HVl to observe 


such a system, a probing field is applied to the principal 
system, whose output is observed, and then we can construct a 
quantum filter. A possible experimental realization is discussed 
in Section |Vl Conclusions are drawn in Section IVll 

II. Dynamics of the ancillary system driven by 

WHITE NOISE 

A. Review of Markovian quantum systems 

1) Hamiltonian: Consider a quantum system defined on 
a Hilbert space f) interacting with an electromagnetic field 
defined on the Boson Fock space ^ over L^(R+). For example, 
this system may be represented by an optical mode in a cavity 
interacting with a probing field. The unitary dynamics of the 
total system on the space 1) 0 5^ is defined by the Hamiltonian 

H = Hs + Hi + Hf (1) 

with the total system evolution operator At = exp(—iiTi), 
where Hs is the system Hamiltonian describing its free 
evolution on the space f) and hereafter we set Planck’s constant 
h — 1. The free evolution of the field on the space ^ is 
described by a unitary operator 0 f = exp(—iiFpf), which is 
determined by a field Hamiltonian Hp = {(jj)b{u!)duj 

with boson annihilation (creation) operator b{uj) (b^{uj)) de¬ 
fined on the space ^ satisfying a delta commutation relation 
[ 6 (a;), = S{uj — ui'). The interaction Hamiltonian Hi 

can be expressed as Hj = i {b'^{w) L — L'^ b{uj))duj, where 

the coupling operator L acting on the space t) only is expressed 
as a product between a decoherence rate ^7 S R+ and a 
system operator. 

2) Dynamical Equation: In the interaction picture, the 
effective Hamiltonian can be obtained as 

iFeff(t) = elHOt = Hs + iib^{t)L - Lh{t)), (2) 

where the field is defined as 

b{t) = —= / (3) 

V J —QQ 

satisfying the delta commutation relations 

\b{t),b^{t')] = 5{t-t'),\b{t)Mt')]=Q- (4) 

The field may be treated as a quantum stochastic process. 
Note that we assume the initial state of the field on the Fock 
space is a vacuum state such that this process is analogous 
to Gaussian white noise with zero mean. 

With the definition of the field Q, we define an integrated 
operator by Bt = J*^b{t')dt', bI = {t')dt' where 

[Bt,Bl,] = mm{t,t'),[Bt, Bt>] = 0 and thus the operator 
Qt = Bt + b\ is the quantum analog of the Wiener process 
and q{t) — b{t) -f b^{t) is quantum white noise. Hence, the 
evolution operator of the total system in the interaction picture 
Ut = OJAj satisfies a quantum stochastic differential equation 
as follows 

dUt = {- {iHs + \L^L)dt + dBlL - L^dBt}Ut (5) 
in Ito form. 


By using Q, a dynamical equation for an arbitrary operator 
X of interest, namely a quantum Langevin equation, can be 
written down as 

dXt = {-i[Xt,Hsit)]+CL,iXt))dt 

+dBl[Xt,Lt] + [LlXt]dBt (6) 

with a generator 

giX) = -i[X,Hs]+CLiX), (7) 

where Xt = UlxUt, Lt = U}LUt, and Hgit) = U^HgUf 
The notation £.{■) defines a Lindblad superoperator which can 
be calculated as Cn{0) = dN^O,N] + ^[NfO]N for two 
arbitrary operators N and O with suitable dimensions. Such 
an equation describes the dynamics of the system driven by 
an external white noise field, which has been widely used in 
the analysis and control of Markovian quantum systems |j9l. 

3) Input-output relations: To observe the dynamics of the 
system, one may consider an output field. The output field 
is the field after interaction with the system, which can be 
defined as i3out(i) = U^BtUt satisfying a QSDE 

d^out(f) = Ltdt + dBt, (8) 

which shows the output field not only carries information of 
the system but also is affected by noise. As a result, the 
output field can be utilized by a quantum filter or a feedback 
controller fT), ll22l . Il2^ . 

4) {S,L,H) description: In considering the interconnec¬ 
tion of Markovian subsystems, an (S', L, H) description has 
been developed. A Markovian system G can be systematically 
denoted as 

G = {S,L,H), (9) 

where the component S is a scattering matrix describing the 
input-output relation of a field passing through beam splitters, 
the operator vector L is a collection of system operators 
interacting with the external fields, and H is the system 
Hamiltonian ll2^ . 

The {S,L,H) description can also concisely describe the 
interconnection among subsystems.When we consider the con¬ 
nection between two subsystems, Gi and G 2 , they can be indi¬ 
rectly connected in a series product Gi<\G 2 or a concatenation 
product Gi ffl G 2 way via the input-output fields ll24l . For a 
class of subsystems, they can also directly interact with each 
other which is denoted as Gi ixi G 2 ll23ll . With these basic 
interconnections, subsystems can be assembled into quantum 
feedback networks UtM or photonic networks ||25]| . 

B. Ancillary system driven by white noise 

To convert white noise to colored noise with a Lorentzian 
spectrum, we consider a Markovian linear quantum system, 
namely an ancillary system, whose mode can be taken as 
an internal mode of the non-Markovian environment. The 
ancillary system is described by Ga = (I, wpajao); 

i.e., an optical mode in a cavity, where wq is the angular 
frequency and oq (aj) is the annihilation (creation) operator 
of the ancillary system defined on a Hilbert space ho- Here 



the coupling operator is chosen as y^oo, where is a 
damping rate with respect to the white noise field. 

Let Ut denote the unitary evolution for the ancillary system 
in the interaction picture with respect to the white noise field, 
which satisfies a QSDE as follows 



+ - y/^aldBt}Ut. {10) 

The generator for the ancillary system is Qa(^a) = 
-i[Xa, wpajao] + £^ao(^a), where is an operator of 
the ancillary system. Note that since the ancillary system is 
driven by white noise, we continue to use the notation Ut and 
dBt to describe the system evolution operator and white noise 
process, respectively. Hence, a QSDE for the operator ag can 
be obtained as 

dao(t) =-(y+ iuJo)ao(t)dt - y/^dBt (11) 

with ao(t) = UtUoUt. 

We define 

c(t) = ( 12 ) 

as a fictitious output, which is different from the output field 
discussed in III- A3 1 Then we have that the fictitious output c(t) 
satisfies a QSDE as follows 

dc(t) = -(y + iuJo)c(t)dt + yds* (13) 

whose formal solution can be expressed as 



Fig. 1. Schematic diagram for the direct coupling between the ancillary and 
the principal system. 

with the memory kernel function 

M(t-t')=J-^lS(u;)]= f at - T)Cit' - T)dT, (18) 

which is different from the delta commutation relation for 
the white noise (IHl. Here, denotes the inverse Eourier 
transform. 

This model has a finite bandwidth which is determined by 
7 o. In the broadband limit; i.e., 70 — 00 , we have 

/ t 

b{t) (19) 

-00 ^ 

where the fictitious output c{t) reduces to white noise with 
delta correlation function. This means the dynamics of the 
ancillary system can be ignored in the broadband limit. 


c{t) = + 

with an initial state c(<o)- 


Jto 2 


(14) 


C. Lorentzian spectrum and broadband limit 

As a part of the environment, the dynamics of the ancillary 
system may be assumed to start from a long time ago in which 
case we let tg —— 00 . Hence we have 

cit) = U^cUt= f ^(f-r)6(r)dT (15) 

J —00 


III. Principal system interacting with the 

ANCILLARY SYSTEM 
A. Dynamics of the augmented system 

A principal system Gp on a Hilbert space 1) with a free 
Hamiltonian Hs is of interest, which is disturbed by the 
colored noise created by the noise model, i.e.,the ancillary 
system, via a direct interaction. Thus the principal system 
and the ancillary system constitute an augmented system. 
We assume that the interaction Hamiltonian for the coupling 
between the principal system and the ancillary system is 

Hi = i{az - Z^c), Z = y/^K, (20) 


in place of (fl4l i. which is a convolution involving the white 
noise field and the kernel ^(f) = The power 

spectral density for c{t) is Lorentzian calculated to be 


S{uj) 



utg)^ 


(16) 


with a center frequency ujg and a linewidth 70 determined by 
the angular frequency of the ancillary system and the damping 
rate with respect to the white noise field, respectively. 

In this case, the commutation relation for c(f) is calculated 
to be 

(17) 


where Z is an direct coupling operator of the principal system 
expressed as a product between a principal system operator 
K and a coupling strength i/k. The principal and ancillary 
systems influence each other due to their direct interaction as 
shown in Eig. [T] This augmented principal-ancillary system 
can be described by using an (5, L, H) description as 

Gp^a = {l,y/^ag,Hs + Hi + UJga\ag), ( 21 ) 

where the evolution operator Ut of the total system satisfies a 
QSDE as follows 

dC7t = { - i{HsH i + ujgalag)dt - Yt^g^odt + 

y/y^dBlag - y/y^aldBt}Uf ( 22 ) 












Let X' denote any operator for the augmented principal and 
ancillary system. Its evolution can be defined as X'^ = u}x'Ut 
which satisfies a QSDE written as 

dX' = + 

+ {[K4Zt] + [Zl-cun)^t 
+y^(d4[X;,do(f)] + [al{t),X[\ABt), (23) 

with Jit = uJHs +UQalao)Ut, do(f) = UjoUt, Ct = 
UjcUt, and = fjIZUt. 

In particular, for X' = X a principal system operator, 
Eq. (l23b reduces to 

dXt = -i[Xt,Hsit)]dt + (4[Xt, Zt] + [Z},Xt]ct)dt, (24) 

with Xt = UjXUt and Hs{t) = UlHsUt- When X' = c = 
— i.e., for an operator of the ancillary system, we have 

dc(f) = -(y + ia;o)c(f)df + ^^tdf + (25) 

Note that c{t) can be written as 

If* 

c(t) = c{t) + - T)ZrdT, (26) 

^ -Ito 

which shows the ancillary system not only depends on the 
Lorentzian noise c{t) defined in (fTdl i but also is disturbed by 
the principal system as indicated by the integral term in (l26l l. 

B. Interaction picture with respect to the ancillary system 

We can move to an interaction picture with respect to the 
ancillary system by defining an evolution operator as V) = 
UjUt, whose evolution satisfies 

Vt = {-iHs- (Z^cit) - c^(t)Z)}Vt. (27) 

In this interaction picture, the system is described by 

G={-,-,Hs + i{cHt)Z-Z^cm (28) 

where c(<) is given by (fTTT i with a Lorentzian spectrum. 

Note that in the interaction picture the system is driven by 
Lorentzian noise as given in Eas. (l27l l and (l28l l. The evolution 
of an operator X for the principal system in the interaction 
picture is equivalent to that in the augmented system due to 

Vt^XVt = UlUtXUlUt = ulXUtUlUt = UlxUf (29) 

Hence, the operator evolution for the principal system in 
Eq. (l24l l is disturbed by Lorentzian noise as well. 

C. Non-Markovian dynamics of the principal system and its 
Markovian limit 

Substituting the solution (l26l l into (l24l i. a non-Markovian 
Langevin equation for the principal system can be obtained as 

Xt = -i[Xt,Hs{t)]+c\t)[Xt,Zt] + [ZlXt]c{t) 

+D{f\Zj[XuZt] + [zlXt]D{^,Z)t (30) 
where the convolution terms are expressed as 

m,Z)t = 11 at-r)Z,dT. (31) 



Fig. 2. Schematic diagram for probing a non-Markovian quantum system. 

This Langevin equation coincides with the existing non- 
Markovian Langevin equations whose integral terms represent 
the memory effect ifTTI . 12^ . 

Remark: Note that taking the broadband limit 70 — 
+c», i.e., using the quantum Wong-Zakai theorem lIZTl . the 
Langevin equation (l30l l reduces to 

X* = -i[Xt,Hs{t)]+C^i,^{Xt) 

+v^64f)[Xt,Xt] + ^[KlXt]h{f), (32) 

with Kt = U^KUt, which is coincident with the Markovian 
Langevin equation 0 . 

D. Master equation 

By using the fact that the expectation of an operator X 
in the Heisenberg picture is equal to that in the Schrodinger 
picture, we can obtain an unconditional master equation for 
the augmented principal and ancillary system as 

pt = -i[Hs+u}oalao,pt]+ 

+ JZ,pt] + [pt,Z^c], (33) 

where pt is the unconditional state of the augmented system 
and the superoperator £*(•) is the adjoint of the Lindblad su¬ 
peroperator calculated as C*^{0) = \N[0, X'l'] + i[X, OjX't' 
for arbitrary operators N and O with suitable dimensions. 

As shown in Eq. (l33l l. the state evolution of the augmented 
principal and ancillary system is Markovian, where the state 
variation only depends on the present state. One can also 
obtain the unconditional state pf^ of the principal system by 
calculating 

Pt=tra[pi], (34) 

which will not satisfy a Markovian evolution. Note that tro[-] 
means the partial trace with respect to the ancillary system. 

IV. Quantum filtering for non-Markovian 

QUANTUM SYSTEMS 

A. The augmented system under a probing field 

To estimate the dynamics of the non-Markovian system, we 
can construct a quantum filter using a probing field defined 











on a Fock space as shown in Fig. |2] The total system Gt 
can be described as 

Gt = (I, ^,Hs + Hi+ uj^ala^) (35) 

where L is the coupling operator of the principal system for 
the probing field. We denote the evolution operator of the total 
system as JJt which satisfies a QSDE as follows 

dUt = { - i{Hs + Hi + uJoalao)dt - ^alaodt 

-^L'^Ldt + ,/^dBlao - 

+dBlL- L'<dBt}Ut. (36) 

Then a QSDE for an operator X' of the augmented principal 
and ancillary systems dehned on f) ® f)o can be derived as 

dl; = -i^, Ht]dt + ([!', dlzt] + [zjdt, xi])dt 

+ ('^V5o5o(t) (-^t) + {X^))dt 

+^idB}[Xi, do(i)] + [4(0: K]dBt) 
+dBl[XiLt] + [Llxi]dBt, (37) 

where 1' = UjX'Ut, Ht = uKHs + tcoajao)!/*, c* = 
illcUt, Zt = UlZUt, It = ill Lilt, do(t) = BjaoBt and 
dBt is the probing held process. 

Note that supposing an operator of the augmented system 
can be denoted as X' = Xp0Xa, the generator can be written 
as 

St(X') = gp(Xp)0Xa+Xp0ga(Xa) 

-ilx', Hi], (38) 

where gp{Xp) = -i[Xp,Hs] + Cl{Xp) and ga{Xa) = 
-i[Xa,Woajoo] + C^aoiXa) are the generators for the 
principal system and the ancillary system, respectively. 

B. Belavkin quantum filter 

Using the probing held, the system can be continuously 
monitored via homodyne detection, where a quadrature of the 
probing held is detected and can be used as input to a quantum 
hlter. 

The output Yt satisfying the QSDE dYt = {Lt+L\)dt+dQt 
with Qt — Bt + B}, commutes with an observable of the 
augmented system, i.e., a non-demolition condition 

[X^,Yr] = 0, 0<T<t, (39) 

is satished which means the continuous measurement of the 
held does not change the observable of the system. Hence, a 
quantum hlter using the output process {Yr,0 < r < t} for 
the estimation of the evolution of the system observable X'^ 
can be constructed based on. Note that we have assumed that 
the detection efficiency of the homodyne detector is perfect, 
i.e., 100% efficiency, and the state of the held is in a vacuum 
state. 

The estimate of an observable X^ is dehned by a conditional 
expectation as 

x', = ^t{x')=nx't\ytl (40) 


where yt is a commutative subspace of operators generated 
by the measurement results Y{t), 0 < t < t. The conditional 
expectation can be interpreted as the orthogonal projection of 
X'f. onto a subspace of yt, which means that X'^ — TTt{X') is 
orthogonal to the measurement subspace yt, i.e., 

E[( 7 rt(X')-X')Ct ]=0 (41) 

for an arbitrary operator Ct in 3 ^* HI, 1281 . 

Hence, we can obtain a Belavkin quantum hlter for the 
augmented system as 

dTTtiX') = 7Tt{gTiX'))dt-{TTtiX'L + L^X')-TTtiX') 
XTrt{L + L^))idYt - TTt{L + L^)dt) (42) 

where dW = dYt — TTt{L + L^)dt and W is called the 
innovation process and is equivalent to a classical Wiener 
process. Note that the increment dW is independent of 
7rT-(X'), 0 < T <t. 

C. Stochastic master equation 

The conditional expectation irfiX') is dehned for the aug¬ 
mented system and thus a conditional density matrix pt for 
the augmented system can be dehned by 

TTt{.X')=tY[ptX']. (43) 

Hence, a stochastic master equation for the augmented system 

can be obtained from the quantum hlter (l42li as 

dpt = g^{pt)dt +FL{pt)dW (44) 

with = Lpt + PtL'^ - tr[(L -f L^)pt]pt, which is a 

Markovian stochastic master equation. The superoperator g^ 
is the adjoint of (/t- 

However, a stochastic master equation for the density opera¬ 
tor pf of the principal system of interest is not in a Markovian 
form. Instead, we can trace out the ancillary system to obtain 
the conditional state of the principal system pf as 

pf = tra[pt]. (45) 

In practice, one cannot obtain an exact description of pf due to 
the inhnite dimensional nature of the ancillary system. While 
a truncation can be made for the ancillary system, i.e, we 
can assume it is a iV-level system and thus it is possible to 
calculate an approximation to the partial trace (l45l l. 

V. An illustrative example 
In this section, an optical realization of a non-Markovian 
quantum system is shown in Eig. [3 The ancillary cavity, which 
is vertically oriented, is pumped by white noise. Another cavity 
is orthogonally oriented to the ancillary cavity and can be 
considered as the principal system. The optical modes in the 
two cavities are directly and strongly coupled by an optical 
crystal. A probing field is applied to the principal cavity, whose 
output is observed via homodyne detection in order to apply 
a quantum filter to the system. 

The Hamiltonian of the principal system can be written as 
Hs = Wsojas with an angular frequency Wg and an annihila¬ 
tion (creation) operator Og (aj). Then, the coupling operator 
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Fig. 3. An illustrative example of optical systems. 


Z and L can be specified as Z and L = 

Substituting these operators into Eq. (iJTl l. Langevin equations 
for the augmented principal and ancillary cavities can be 
written as 

hs(t) 

_ ao{t) 


where bp{t) and b{t) are the probing field and white noise field, 
respectively. The output equation with respect to the probing 
field bp{t) is 


, 7i 

2 2 


ttsit) 

v'kTo . 70 

2 ^^0 2 


_ ao{t) 


V^bpit) 

V^b{t) 


,(46) 


bout{t) = bp{t) + y/^as{t). (47) 

Since the operators in Eq. (l46l l are not self-adjoint operators 
and the coefficients may be complex valued, it is convenient 
to express Eq. (l46l l and WT\ as 


x{t) = Ax{t) -b Bu{t) 

y{t) = C [ / 0 ] x{t) + [ r 0 ] u{t) 

with 


71 

2 

UJs 2 

0 

-UJs 

_2l 0 

v'kTo 

2 

2 

0 

^ 2 

uJo 

0 

\/^ 

2 

70 

2 -1 

diag[-v7i, 

-\/7i,-\/7o, - 

-vTo], 


C = diag[y^, 


(48) 

(49) 


(50) 

(51) 

(52) 


in a quadrature representation with self adjoint 
operators and real valued coefficients, where 
x{t) = [qs{t),Ps{t),qo{t),po{ty\^, u(t) = 

[vp{t),Vg{t),Wp{t),Wg{t)]'^, and y{t) = [yp{t),yq{t)]'^ are 
the quadrature representations of the operators of the principal 
and ancillary systems, the probing and white noise fields, 
and the outputs of the probing field, respectively. By using a 

" 1 1 


transformation matrix ^ 

v2 

x{t), u{t), and 


^ ^ , the components of 

y{t) are calculated as [qs(t),Ps(t)]’^ = 

[qoit),Po{t)f = 5[ao(f),a;(f)]'^, 


[Vp{t),Vq(t)f = ^[bp{t),bl{t)f, [Wp{t),Wg{t)f = 

E[b{t),b^{t)f,[yp{t),yg{t)f = 

A. The mean of the principal system 

Assume the unconditional state of the augmented system is 
Gaussian ll29l and thus the means of operators for the system 
(|48] | can be used to observe the system dynamics. Hence, we 
consider m{t) = {x{t)) satisfying 

m(f) = Am(t), (53) 


where the quantum expectation (•) = tr[-po] is with respect to 
the initial state of the augmented system pQ. 

On the other hand, the quantum filter (l42l l as the estimation 
of the principal and ancillary systems is actually a quantum 
Kalman filter due to the linearity of the augmented sys¬ 
tem ]3oi, im. The conditional system dynamics are governed 
by the quantum Kalman filter as follows 


dxt 

% 


with E = 


© 


Axtdt + {VtF^ + (54) 

AVt + VtA^ + D- (VtF^ + S^Im(Ff) 
x(FVt+Im(E)S) (55) 

=1 ° n , L = Fx=lyW^-\0]x, F = 


F+F*, and D = T,Re{E^F)T,'^, where x(t) is the conditional 
expectation of x{t) and V) is a symmetrized covariance matrix. 
The bold 0 is a 2 x 2 zero matrix. 

Eor simplicity, we choose the parameters of the system as 
uJs = ujQ = lOGHz, K = 0.6, 7 o = 0.6, and 71 = 0.8 
and assume that the initial mean m( 0 ) = [ 1 , 0 , 0 , 0 ]^ of the 
unconditional state is the same as the initial conditional state 
expectation of the Kalman filter, i.e., x( 0 ) = [ 1 , 0 , 0 , 0 ]^ where 
the first element of x(0) is (gs(0)) = 1. In Eig. 01 the red 
solid line is the trajectory of the mean of the unconditional 
state qs{t) for the principal system. The oscillations of the 
curve envelopes are caused by the disturbance of the ancillary 
system, which indicates that energy is exchanged between the 
principal and the ancillary system showing non-Markovian 
characteristics. Compared with the unconditional trajectory, 
the average trajectory of the conditional expectation of position 
for 1000 realizations denoted as {qs{t)), is plotted as the blue 
dashed line in Eig. 0] This shows that the blue line for the 
quantum filter can match with the red line for the unconditional 
state of the system, i.e., the quantum filter can estimate the 
unconditional state of the system. 


B. The spectrum with respect to the white noise field 

Alternatively, it is useful to observe the output field spec¬ 
trum as considering the properties of the system EH. To 
calculate the power spectral density of the output of the prob¬ 
ing field, we can transform (06l l and (07l i into the frequency 
domain via the Eourier transform and express the output in 
terms of the white noise field and the probing field as 

(>out(iw) = Gi{\uj)bp{\uj) F G 2 (ia;) 6 (ia;) (56) 

in the frequency domain, where 6 out(iw), bp{iuj) and 6 (ia;) are 
the Eourier transform of 6 out(f), bp{t) and b{t), respectively. 







































Fig. 4. The dynamics of the unconditional and conditional means of the 
position components for the principal system with k, = 0.6, A = 0,70 = 0.6, 
and 7i = 0.8. 


Fig. 5. The spectrum \G 2 {<^)\^ varies with the coupling strength k, where 
A = 0,70 = 0.1, and 71 = 0.8. 


The transfer function with respect to the probing field and the 
white noise field can be calculated as 

(i(a;s - g;) - 7)(i(^o " ^) + t) + ^">'0 


Gi(iw) = 


G2(iw) = 


(i(ws - w) + 7)(i(‘^o - w) + 7) + ^^70 ’ 

_-7o7^_ 

(i(ws - w) + 7)(i(wo - w) + 7) + «7o ’ 


where Hj = ojg — OJ, A = lOs — ujq, and 

T(w) + 77o-^)" 


IGiMp 

\G2{A)\^ 

_ /7i^ 


T(w) + (k7o ■ 


7071 12 
4 / 


«^7i7o 


T((7) + 770 + 277)2 


(57) 


(58) 


respectively. 

By detecting the quadrature of the output field 7fout(iaj) = 
5 [feout(iw) + 6out(i‘7’)], the power spectral density of the output 
field can be calculated as 

SiCo) := (|Vo„t7)p) = ' 



(|Gi7)p + |G 2 (t:;)n (59) 


(60) 


with T(ai) = (2i- + — A)^ + AA- 2K^oA{oj — A). 

Although the total output spectrum is flat as given in (l59T l 
due to the passivity properties of the system |[32l, we can apply 
a coherent probing field whose strength is much higher than 
the quantum white noise and thus the spectrum |Gi(w)p can 
be observed so as to calculate the spectrum |G 2 (£l))P which 
reflects the influence of the ancillary system on the output 
held. 

Fig.E] shows the power spectral density |G 2 ((I')P varying 
with the coupling strength k. Here, we assume there is no 
detuning A = 0 and the principal system damping rate 
7 i = 0.8 is much higher than that with respect to the white 
noise 70 = 0.1. As the coupling strength k is increased, the 
amplitude of the noise at the system working frequency is 
decreased while the bandwidth of the spectrum is broader, 
which means the coupling strength k can affect both the noise 
amplitude and the bandwidth. 

The power spectral density |G 2 (w)p varying with the detun¬ 
ing A is plotted in Fig. | 6 ] with parameters k = 0.1 ,70 = 0.1, 


Fig. 6. The shape of 102(1^)12 varies with A, where re = 0.1,70 = 0.1, 
and 71 = 0.8. 


and 7 i = 0.8. When there is no detuning as given by the 
blue line, the noise is strong at the system frequency. As the 
detuning is increased via decreasing the angular frequency of 
the ancillary system, the spectrum |G 2 (a’)p is driven away 
from the system frequency, whose amplitude is decreased as 
well. This illustrates that the non-Markovian effect generated 
by the ancillary system becomes weaker as the detuning is 
increased. 

Fig. |7] shows the power spectral density |G 2 (A)p can also 
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Fig. 7. The spechum |G2(ta)|2 for vaiying 70 with A = 0, re = 0.1, and 
71 = 0.8. 













































be varied with the damping rate 70 with respect to the white 
noise. As predicted, the bandwidth of the Lorentzian spectrum 
is broader as 70 increases, varying spectrum |G' 2 (w)P with 
respect to 70 gives the same result, which indicates that 70 
determines the bandwidth of the noise spectrum. 

VI. Conclusion 

In this paper a Markovian representation approach to mod¬ 
elling a non-Markovian quantum systems compatible with 
quantum filtering theory has been presented. The ancillary 
system of this model plays the role of the internal mode of 
the environment, which can convert white noise to Lorentzian 
noise, resulting in non-Markovian dynamics of the principal 
system. In addition, the quantum filter has been derived for the 
non-Markovian quantum system. Simulation results show the 
quantum filter can estimate the state of the non-Markovian 
quantum system. Since the total system is expressed in an 
extended Hilbert space, this is equivalent to using a Markovian 
network to describe both the non-Markovian system and its 
environment. For future work, such an approach could be 
extended to represent a non-Markovian quantum system with 
arbitrary colored noise by designing the structure of the 
Markovian network such that a more general quantum filter 
can be constructed for the non-Markovian quantum systems. 
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